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We present branching fraction measurements for the decays B° — > p°K*° , B° — > foK*° , and B° — >• 
p~K* + , where K* is an S-wave (Ktv)q or a K*(892) meson; we also measure B° — > foK% (1430)°. 
For the K*(892) channels, we report measurements of longitudinal polarization fractions (for p 
final states) and direct CP-violation asymmetries. These results are obtained from a sample of 
(471.0 ± 2.8) x 10 6 BB pairs collected with the BABAR detector at the PEP-II asymmetric-energy 
e + e" collider at the SLAC National Accelerator Laboratory. We observe p°-fC(892)°, p°{K-k)* q °, 
/o7C(892)°, and p~A'*(892) + with greater than 5a significance, including systematics. We report 
first evidence for fo(Kir)Q° and /oA'|(1430) , and place an upper limit on p~ (Ktv)q + . Our results 
in the if* (892) channels are consistent with no direct CP violation. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



I. INTRODUCTION 

Measurements of the branching fractions and angular 
distributions of B meson decays to hadronic final states 
without a charm quark probe the dynamics of both the 
weak and strong interactions. Such studies also play 
an important role in understanding CP violation in the 
quark sector and in searching for evidence for physics 
beyond the Standard Model 

We report measurements of branching fractions for 
the decays B° -> p°K*{892)°, B° ->■ f a K*(892)°, 
B° -> p-K*(892)+, B° -> p°(Kny °, B° -> / (iv7r)*°, 
B° -> p-(Kir)* + , and B° -> f K*(U30)°. For the 
pK* (892) channels we measure the longitudinal fraction 
/l, and for all if* (892) channels we measure charge 
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asymmetries A c h- The notation p refers to the p(770) 
and /o to the / (980) 3. Throughout this paper we 
use K* to refer to any of the scalar (Ktt)q, vector 
if* (892), or tensor iC*(1430) states 0. The notation 
(Ktt)q refers to the scalar Kw, which we describe with 
a LASS model d i], combining the ^(1430) reso- 
nance with an effective-range non-resonant component. 
Charge-conjugate modes are implied throughout this pa- 
per. 

The charmless decays B — ¥ pK* proceed through dom- 
inant penguin loops and CKM-suppressed tree processes 
(B + —> p + K*° is pure penguin), as shown in Fig. [T] 
Naive factorization models predict a large longitudinal 
polarization fraction fi, (of order (1 — Arn^/m 2 B ) ~ 0.9) 
for vector-vector (VV) decays, where my and tub are 
the masses of the vector and B mesons, respectively 
However, measurements of penguin-dominated VV de- 
cays, such as the previous measurements of B — > 
p°K* (892)° and B+ -> p+K* (892)° 0,0, find f L - 0.5; 
a recent BABAR measurement of B + p°K *(892)+ finds 
f L = 0.78 ± 0.12 (p. Recent predictions in QCD Factor- 
ization (QCDF) .9] can accommodate Jl ~ 0.5, although 
correctly predicting both the branching fraction and /l 
remains a challenge. 

Both the BABAR. and Belle Collaborations have pre- 
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the notation % r 



cos(# p o). 



(d) 



FIG. 1: Feynman diagrams for (a-b) B° — > p°K*° and (c-d) 
B° ->■ p-K*+. Gluonic penguin diagrams (a, c) dominate 
over tree (b, d) contributions. 



viously measured the branching fractions of B° — > 
p° K* (892)° and P° -> / if*(892)°. BAR4R has also 
placed a 90% confidence level (C.L.) upper limit on 
B° -> p- if* (892)+ @, Q. Belle searched for non- 
resonant P° — > p°K + ir~ and JqK + tt~ decays, find- 
ing a five standard deviation (5er) significant result for 
p°K + ir~ 0]. Decays involving a p or /o along with a 
(^7t)q or K% (1430) have not been the subject of previous 
studies. Predictions exist from both QCDF and per- 
turbative QCD (pQCD) [HI for the branching fractions 
(B) of the pX *(1430) channels, with QCDF predicting 
values C(few x lO" 5 ) and pQCD 0(5 x 10~ 7 - 10~ 5 ). 
Improved experimental measurements will help refine 
predictions and constrain physics beyond the Standard 
Model. 

The decays B° -> p°K* (892)° and B° -> p~K* (892)+ 
are of the form B — > V V; these decays have three polar- 
ization states, which are, in principle, accessible experi- 
mentally. In practice, a full angular analysis requires a 
large number of signal events. In the analyses described 
in this paper, we integrate over the azimuthal angle (the 
angle between the two vector meson decay planes). The 
azimuthal angle is not correlated with any specific direc- 
tion in the detector, so we assume a uniform acceptance 
over this angle. We define the helicity angles 9 k* and 
9 P and the azimuthal angle cf> as shown in Fig. [5] The 
helicity angles are defined in the rest frame of the vector 
meson: 9k* is the angle between the charged kaon and 
the B meson in the K* rest frame; 9 P is the angle between 
the positively charged (or only charged) pion and the B 
meson in the p rest frame. In the analysis of the f£T*(892) 
channels, we make use of the helicity observables, defined 
for a = p, K* as T~L a = cos(9 a ). Occasionally, we refer 
to a specific charge state, e.g. p° , which we indicate with 




FIG. 2: Definition of the helicity angles for B° -> p°K*°. 

The longitudinal polarization fraction for B° — > 
pK* (892) can be extracted from the differential decay 
rate, parameterized as a function of Ok* and 9 p : 



d 2 T 



r d cos 9 K*d cos 9 p 

j(l - Jx) sin 2 9 K * sin 2 9 P + f L cos 2 9 K * cos 2 9 P 
The CP-violating asymmetry is defined as 

r- - r+ 



(1) 



r- + r+ 



(2) 



where the superscript on the decay width T refers to the 
charge of the kaon from the K* decay. 

All results in this paper are based on extended max- 
imum likelihood (ML) fits as described in Section IVII 
In each analysis, loose criteria are used to select events 
likely to contain the desired signal B decay (Sec. IIIII 
[V]) . A fit to kinematic and topological discriminating 
variables is used to differentiate between signal and 
background events and to determine signal event yields, 
CP-violating asymmetries, and longitudinal polarization 
fractions, where appropriate. In all of the decays ana- 
lyzed, the background is dominated by random particle 
combinations in continuum (e + e _ — > qq, q = u,d,s,c ) 
events. Although qq background dominates the selected 
data sample, background from other BB decays tends to 
have more signal-like distributions in the discriminating 
variables. The dominant BB backgrounds are accounted 
for separately in the ML fit, as discussed in Sec. lIVrJl Sig- 
nal event yields are converted into branching fractions via 
selection efficiencies determined from Monte Carlo (MC) 
simulations of the signal as well as auxiliary studies of 
the data. 



II. DETECTOR AND DATA 

For this analysis we use the full BABAR dataset, col- 
lected at the PEP-II asymmetric-energy e + e~ collider 
located at the SLAC National Accelerator Laboratory. 
The dataset consists of (471.0 ±2.8) x 10 6 BB pairs orig- 
inating from the decay of the T(AS) resonance, produced 
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at a center-of-mass (CM) energy y/s — 10.58 GeV. This 
effectively doubles the dataset from the previous BABAR 
measurement [f|. 

The asymmetric beam configuration in the laboratory 
frame provides a boost to the T(AS) of (3j = 0.56. 
This results in a charged particle laboratory momentum 
spectrum from B decays with an endpoint near 4 GeV. 
Charged particles are detected and their momenta mea- 
sured by the combination of a silicon vertex tracker, con- 
sisting of five layers of double-sided detectors, and a 40- 
layer central drift chamber, both operating in the 1.5-T 
magnetic field of a solenoid. For charged particles within 
the detector acceptance, the average detection efficiency 
is in excess of 96% per particle. 

Photons are detected and their energies measured by 
a CsI(Tl) electromagnetic calorimeter (EMC). The mea- 
sured 7T° mass resolution for 7r°'s with laboratory mo- 
mentum in excess of 1 GeV is approximately 8 MeV. 

Charged particle identification (PID) is provided by 
the average energy loss in the tracking devices and by 
an internally reflecting ring-imaging Cherenkov detector 
(DIRC) covering the central region. Additional infor- 
mation that we use to identify and reject electrons and 
muons is provided by the EMC and the detectors in- 
stalled in a segmented solenoid flux return (IFR) . The 
BABAR detector is described in detail in Ref. [ll[ . 

III. CANDIDATE RECONSTRUCTION AND B 
MESON SELECTION 

We reconstruct B-daughter candidates through their 
decays p° -> 7r + 7r~, / -> 7r+7r~ , p~ ->■ ir~ir°, K*° -> 
K + n~ , K* + — > K + ir°, and ir° — > 77. We apply the same 
selection criteria for /o and p° candidates. 

The if* (892) channels are analyzed separately from 
the (Ktt)q and i\2(1430) decays, though the analyses 
share many similarities, including most event selection 
requirements. Where the analyses differ, we specify 
the If* (892) channels as the "low mass region" (LMR), 
distinguished by the Kn mass requirement of 750 < 
m K -K < 1000 MeV. The {Ktt)* and TV? (1430) anal- 
yses are performed in the "high mass region" (HMR), 
1000 < m & < 1550 MeV. 

The invariant masses of the B-daughter candidates 
must satisfy the following requirements: 120 < rn 77 < 
150 MeV, and either 470 < m T7r < 1070 MeV (LMR) or 
470 < m W7r < 1200 MeV (HMR). The tttt and Kn mass 
intervals are chosen to include sidebands large enough to 
parameterize the backgrounds. 

All photons are required to appear as a single clus- 
ter of energy in the EMC, not matched with any track, 
and to have a maximum lateral moment of 0.8. We re- 
quire the energy of the photons to be greater than 50 
MeV and the ir° energy to be greater than 250 MeV, 
both in the laboratory frame. All charged tracks are re- 
quired to originate from within 10 cm of the beamspot 
in the direction along the beam axis and within 1.5 cm 
in the plane perpendicular to that axis. Charged kaon 



candidates are additionally required to have at least 12 
hits in the drift chamber and a transverse momentum 
of pt > 100 MeV. The charged tracks are identified as 
either pions or kaons by measuring the energy loss in 
the tracking devices, the number of photons recorded by 
the DIRC, and the corresponding Cherenkov angle; these 
measurements are combined with information from the 
EMC and the IFR, where appropriate, to reject electrons, 
muons, and protons. 

When reconstructing p~ and K* + candidates, the 
mass of the 7r° candidate is constrained to its nominal 
value [2j. The ir° is constrained to originate from the 
interaction point, taking into account the finite B meson 
flight distance; the charged track is required to originate 
from the interaction point. For p° and K*° candidates, 
the two charged tracks are required to originate from 
a common vertex, as determined by a generalized least 
squares minimization using Lagrange multipliers; we re- 
quire the change in \ 2 between two successive iterations 
in the fitter to be less than 0.005, with a maximum of 
6 iterations. The B meson candidate is formed by per- 
forming a global Kalman fit to the entire decay chain. 

A £?-meson candidate is characterized kinematically by 
the energy-substituted mass ities arL d the energy differ- 
ence AE, defined in the T(4S) frame as 



mES = y 4 s _ Pb and 
AE = E* B -^, 

where q* B = (E B ,p B ) is the four momentum of the B- 
candidate in the T(AS) frame and s is the square of the 
invariant mass of the electron-positron system. m B s and 
AE are favorable observables because they are nearly 
uncorrelated. The small correlation is accounted for in 
the correction of the fit bias (see Sec. IIXj) . Correctly re- 
constructed signal events peak at zero in AE and at 
the B mass [2| in m B s, with a resolution in m B s of 
around 2.5 MeV and in AE of 17-37 MeV. We select 
events with 5.26 < m ES < 5.2893 GeV. For p°K*°, we 
require |A.E| < 0.10 GeV, while for p~K* + , we allow 
-0.17 < AE < 0.10 GeV to account for a long low-side 
tail resulting from poorly reconstructed 7r°'s. 



IV. SOURCES OF BACKGROUND AND 
SUPPRESSION TECHNIQUES 

Production of BB pairs accounts for only about 25% of 
the total hadronic cross section at the T(AS) peak. The 
bulk of the cross section arises from continuum events. 
Tau-pair production and other QED processes contribute 
as well. We describe below the main sources of back- 
ground and discuss techniques for distinguishing them 
from signal. 
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A. QED and tau-pair backgrounds 

Two-photon processes, Bhabha scattering, muon- and 
tau-pair production are characterized by low charged 
track multiplicities. Bhabha and muon-pair events are 
significantly prescaled at the trigger level. We further 
suppress these and other tau and QED processes via a 
minimum requirement on the event track multiplicity. 
We require the event to contain at least one track more 
than the topology of our final state. These selection cri- 
teria are more than 90% efficient when applied to signal. 
From Monte Carlo simulations [l2[ we determine that the 
remaining background from these sources is negligible. 

B. QCD continuum backgrounds 

The dominant background arises from random com- 
binations of particles in continuum e + e~ — > qq events 
(q = u,d,s,c). The angle 6t between the thrust axis 
[13j of the B candidate in the T(4S) rest frame and that 
of the remaining particles in the event is used to sup- 
press this background. Jet-like continuum events peak 
at values of | cos#t| close to 1, while spherical B decays 
exhibit a flat distribution for this variable. We require 
that events satisfy | cos#t| < 0.7. 

Further rejection is achieved by restricting the range 
of the helicity angle T-L of the p and K* mesons (see 
Fig. 02). We require \H p o\ < 0.9, -0.8 < H p+ < 0.9, 
-0.85 < H K *° < 1-0, and -0.8 < H K *+ < 1.0. These 
requirements reject regions of phase space with low mo- 
mentum 7r + 's and 7r°'s, where backgrounds are typically 
large. 

Additional separation of signal and background is pro- 
vided by a Fisher discriminant J- exploiting four variables 
sensitive to the production dynamics and event shape: 
the polar angles (with respect to the beam axis in the 
e + e~ CM frame) of the B candidate momentum and of 
the B thrust axis; and the zeroth and second angular mo- 
ments Lq,2 of the energy flow, excluding the B candidate. 
The moments are defined in the CM frame by 

Lj = ^Pi |cos0i| J , (3) 

i 

where i labels a track or EMC cluster, &i is its angle with 
respect to the B thrust axis, and pi is its momentum. 

We find that T in continuum background is mildly cor- 
related with the tagging category [14j . which identifies 
the flavor of the other B in the event and places it into 
one of six categories based upon how it is identified. Al- 
though the tagging category is not used elsewhere in this 
analysis, we find that the overall signal-to-background 
separation provided by T can be slightly improved by 
removing this correlation. For each tagging category as 
well as the category for which no B tag is assigned, we fit 
the T distribution with a Gaussian with different widths 
above and below the mean. We then shift the mean of 
the J 7 distribution in each tagging category to align it 



with the average value of the T means in all tagging cate- 
gories. The T distributions typically have a mean around 
—0.25 with an average width around 0.45; shifts are less 
than 0.03 for all categories except for the lepton-tagged 
events (the tagging category with the highest purity), 
for which the shift is about 0.35. The Fisher variable 
provides about one standard deviation of discrimination 
between B decay events and continuum background. 

C. B — > charm backgrounds 

We suppress the background from B mesons decay- 
ing to charm by forming the invariant mass from 
combinations of two or three out of the four daughter 
particles' four- momenta. For p°K*° , we consider D can- 
didates decaying to K~ir + and K~n + ir + . For p~K* + , 
we consider the combinations K~ir + and K~Tr + ir°. The 
event is retained only if \mjj — mp DG \ > 40 MeV for all 
cases except for the D meson formed with K + n~ in the 
p~K* + channel, where we require \mo — m^ 00 ] > 20 
MeV; m^) DG is the nominal D + or D° meson mass 0. 

These D vetoes greatly reduce the amount of B — > 
charm background in our samples, but as many of these 
channels have large branching fractions ©(lO^ 1 — 10~ 3 ), 
we include several charm backgrounds as separate com- 
ponents of the maximum likelihood fit, as detailed in 
Sec. lTVDl 

D. BB backgrounds 

Although the dominant background arises from contin- 
uum qq events, care must be taken to describe the back- 
grounds from other B decays, as they have more signal- 
like distributions in many observables. For p°K*(892)°, 
we consider seven BB background categories: B° — > 
p (Kw)* °; B° -> f (Kw)* °; B° -> / 2 (1270)iT(892)° 
with / 2 (1270) -> 7T+7T-; B° -> K+ with -> p Tr-; 
B° —> D~ir + with D~ —> K + ir~ir~\ a combination of 
three B — s- D°X channels with D° -> K+tt'tt ; and 
a branching fraction- weighted combination of 13 other 
dominant charmless B decay channels (charmless cock- 
tail), which have a high probability of passing our selec- 
tion. The dominant channels in the charmless cocktail 
are B + -> a\K + with a" -> p~n + and B + -t rj'K + with 
rf — > p°j. Most channels in the cocktail include a real p° 
or K*°. The number of expected events in each category 
is given in Table HI 

For the {Ktt)^° and / ^ 2 *(1430)° signals, the back- 
ground categories are the same, except that the K* (892)° 
replaces the {Ktt)q° in the first two background cate- 
gories. As will be described in Sec. IVIB1 the first stage 
of the fit to the HMR is insensitive to the composi- 
tion of the 7r + 7r - mass spectrum; therefore p°K *(892)°, 
f K*(892f, and /2(1270)if*(892) are included in the 
same K* (892)° category. Additionally, due to the wider 
K + it~ mass range in the HMR, 28 charmless B decay 
channels are combined in the charmless cocktail. 
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TABLE I: BB background categories for B° -> p°K*° and 
expected yields in the LMR and HMR. 



p°K*° background LMR HMR 

B° -> p°(Ktt)o° 215 ±34 

B° -> fo{K-K)l° 19 ±6 

B° -> / 2 (1270) JT (892)° 47 ± 3 

B° -> aT/if* 15 ±3 40 ± 9 

B°^D"7r + 209 ± 10 922 ±45 

S ->■ _D°X 433 ± 23 1798 ± 83 

charmless cocktail 76 ± 22 149 ± 34 



event in data is 1.02 for p°/f K*° and 1.16 for p~K*+. 
We select the candidate with the highest x 2 probability in 
a geometric fit to a common B decay vertex. In this way 
the probability of selecting the correctly reconstructed 
event is a few percent higher with respect to a random 
selection. 

The sample sizes for the decay chains reported here 
range from 9700 to 37 000 events, where we include side- 
bands in all discriminating variables (except the helici- 
ties) in order to parameterize the backgrounds. 



VI. MAXIMUM LIKELIHOOD FIT 



In analyzing p~K*(892) + , we consider four BB back- 
ground categories: B° — > p~ (Ktt)* q + , B° — > K + with 
Oj -t p~ir°, B a -t p + p~ , and B~ -> D°p~ with 
D° — > K + ir~ir°. The number of expected events in each 
category is given in TableHH For the HMR, p"X*(892)+ 
replaces the signal mode p~(Ktt)q + as a background; the 
other categories remain the same. 



The candidates that satisfy the selection criteria de- 
scribed in Sees. IlllflVl are subjected to an unbinned, ex- 
tended maximum likelihood fit to extract signal yields. 
In all fits, the signal and BB background components 
are modeled with a Monte Carlo simulation of the decay 
process that includes the response of the detector and 
reconstruction chain 1121. 



TABLE II: BB background categories for B° ->■ p~K* + and 
expected yields in the LMR and HMR. 



p-K" 


f background 


LMR 


HMR 


B° — > 


P -(Kw)*+ 


60 ±23 




B° —> 




7± 2 


13 ±3 


B° — !> 


p + p~ 


9± 1 


15 ±2 


B- -> 


D°p~ 


129 ± 17 427 ± 58 



In the HMR fits, the if* (892) yields are allowed to 
float. The HMR (-Ktt)q yields are extrapolated into the 
LMR using a ratio of LMR to HMR MC efficiencies; these 
(Ktt)q background yields are then fixed in the LMR fits. 
The B° -> / 2 (1270) J Ff*(892)° background yield is de- 
termined using a high m T + T - sideband, as described in 
Sec- IVHll thisjield is fixed in the p°K* (892)° fit. 

All other BB backgrounds are modeled from the sim- 
ulation, with yields fixed to experimentally measured 
branching fraction B values 0. For a few channels en- 
tering the charmless cocktail, no B measurements ex- 
ist; in those cases, theory predictions are combined with 
other estimates and a 100% uncertainty is assigned to 
the branching fractions. These unmeasured charmless 
channels account for approximately 26% of the charmless 
cocktail background in the LMR and 40% in the HMR 
(see Table [TJ) . Uncertainties on the BB branching frac- 
tions are accounted for as systematic uncertainties (see 
Sec. ED . 

V. FINAL SAMPLE CRITERIA 

After all selection criteria discussed in Sec. IIIIIITVl have 
been applied, the average number of combinations per 



A. Low mass region fit 

In the low mass region, we obtain the yields, charge 
asymmetries A c h, and longitudinal polarization fractions 
Jl from extended maximum likelihood fits to the seven 
observables: AE, mes, J 7 , and the masses and helicities 
of the two resonance candidates (m^, "niKm Hp, and 
Hk')- The fits distinguish among several categories: qq 
background, BB background (see Sec lIVD]) . and signal. 
The signals p°K* (892)° and / A'*(892)° are fit simulta- 
neously. For each event i and category j we define the 
probability density functions (PDFs) as 

V) = VjirriEs^VjiAE^iF) (4) 
with the resulting likelihood C: 

£ = -^IIE^ (5) 

i=i j 

where Yj is the fitted yield for category j and N is the 
number of events entering the fit. For the p /fo anal- 
ysis, we use the absolute value of Hp in the fit, as the 
distribution is symmetric. We split the yields by the fla- 
vor of the decaying B meson in order to measure A c h- 
We find correlations among the observables to be occa- 
sionally as high as 30% in simulations of the BB back- 
grounds, whereas they are small in the data samples, 
which are dominated by qq background. In signal, cor- 
relations are typically less than 1% and occasionally as 
large as 14%. Correlations amongst observables are ac- 
counted for by evaluating the fit bias (see Sec. IIX|) . 
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B. High mass region fit 

In the high mass region, the ML fit uses the five ob- 
servables: AE, tties, J 7 , m^^, and rriKn- F° r /°~(A7t)q + , 
these five observables are combined in an extended ML 
fit, as above. 

For the p Q /f (Kn)^ and / Af(1430)° channels, we 
perform the ML fit in two stages. Due to the poten- 
tial complexity of the resonant and non-resonant struc- 
tures in the ir + n~ and A + 7r~ invariant mass spectra, 
as well as the fact that many of these structures are 
quite broad, non-trivial correlations exist between sev- 
eral of the ML fit hypotheses. Attempts to perform the 
fit in a single stage using simulated data (see Sec. HXI for 
the general procedure) demonstrate unacceptable con- 
vergence rates in some scenarios. Removing m„ from 
the ML fit greatly improves the convergence rates. We 
therefore employ a two-stage procedure for these HMR 
fits. In the first step, we perform an ML fit using only 
AE, rriES, J 7 , and tuk-k] this allows us to separate out 
"inclusive" (Ktt)q and Af(1430)° signal from qq and 
BB backgrounds. If we observe sufficient (greater than 
3<r statistical significance) signal in the "inclusive" A* 
channels, we perform a second-stage ML fit to m W7r for 
selected signal events. Technical details are given below. 

The PDF for the first-stage fit can be written as 

V) = Vj (mEs^Vi (AE' l )V 3 (r)V 3 (m** J (6) 

for event i and category j . 

In the event of significant signal in the "inclusive" 
(Att)o° or Iff (1430)° channels, we apply the sPlot tech- 
nique [HI to the results of this first fit, which allows us to 
calculate a weight value for each event in each category 
(signal, BB background, etc.) based upon the covariance 
matrix from the likelihood fit and the value of the PDF 
for that event. Specifically, the s Weight for event i of 
category n is given by 



T Nc V 
2^7=1 Vi 



-pi 



(7) 



where N c is the number of categories in the fit, V n j is the 
covariance matrix element for categories n and j, and Y% 
is the yield of category k. 

The s Weight for a given event indicates how much 
that event contributes to the total yield in that category; 
s Weights can be less than zero or greater than one, but 
the sum of all s Weights for a given category reproduces 
the ML fit yield for that category. 

The s Weights from this procedure are used to create 
two datasets: the sWeighted (Att^ and Af (1430)° sig- 
nal samples. These weighted datasets allow us to de- 
termine the 7r + 7r _ mass distribution for the two signal 
samples of interest; these sPlots are faithful representa- 
tions of for the (Ktv)q° and Iff (1430)° signal com- 
ponents, assuming no correlation between and the 
observables used to generate the s Weights. For signal 
MC, we find a maximum correlation of 8% between m vlt 



and the other observables, with correlations typically less 
than 2%. 

In the second stage, we fit the sWeighted (Avt)q and 
K 2 (1430)° m W7r distributions to p° and fo hypotheses. A 
non-resonant it + tt~ component is found to be consistent 
with zero. A er//o(600) component is considered in stud- 
ies of systematic uncertainties (see Sec. lXip . This fit gives 
us the final signal yield for the p°(Ktt)q°, f {K-K)* °, and 
/o-fC!; (1430)° channels. This procedure also determines 
the p° Af (1430)° yield but, as we do not include helicity 
information in the fit, we cannot measure /l, and thus 
we consider that channel a background. 

Due to the two-stage nature of the p° / fo(K tt)q° and 
/oA'j (1430)° fits, the statistical uncertainty has two com- 
ponents. The first is from the uncertainty on the m W7r fit 
to extract the fraction of p°/fo events in the sWeighted 
sample. The second is a fraction of the uncertainty on the 
"inclusive" (Att)* (A 2 *(1430)°) yield, the coefficient of 
which is given by the ratio of p° or fo events to the total 
number of inclusive (Ktt)q° (Af (1430)°) signal events. 



VII. SIGNAL AND BACKGROUND MODEL 

PDF shapes for the signals and BB backgrounds are 
determined from fits to MC samples. For the qq cat- 
egory we use data sidebands, which we obtain by ex- 
cluding the signal region. To parameterize the qq PDFs 
for all observables except ities, we use the sideband de- 
fined by wies < 5.27 GeV; to parameterize m^s, we re- 
quire \AE\ > 0.06 GeV for p°K*° or AE < -0.12 and 
AE > 0.08 GeV for p~K* + . The excluded AE region 
is larger for p~K* + due to the poorer AE resolution re- 
sulting from having two 7r°'s in the final state. 

Signal events selected from the MC contain both cor- 
rectly and incorrectly reconstructed i?-meson candidates; 
the latter are labeled "self-crossfeed" (SXF) . SXF occurs 
either when some particles from the correct parent B me- 
son are incorrectly assigned to intermediate resonances 
or when particles from the rest of the event are used in 
the signal B reconstruction. The fraction of SXF events 
ranges from 2-7% for p° / /oA*° candidates and from 13- 
22% for p~K* + candidates. We include both correctly 
reconstructed and SXF signal MC events in the samples 
used to parameterize the signal PDFs. 

We use a combination of Gaussian, exponential, and 
polynomial functions to parameterize most of the PDFs. 
For the ities distribution of the qq background com- 
ponent, we use a parameterization motivated by phase- 
space arguments [16j. 

In the A* (892) (LMR) fits, the following observables 
are free to vary: the signal yields, longitudinal fraction 
/l for pK* (892), and signal charge asymmetries A c h', 
the qq background yields and background A c h] and the 
parameters that most strongly influence the shape of the 
continuum background (the exponent of the phase-space- 
motivated ui es function; dominant polynomial coeffi- 
cients for AE, resonance masses, and helicities; fraction 
of real p, fo, and A* resonances in the background; and 
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the mean, width, and asymmetry of the main Gaussian 
describing J 7 ). For the HMR fits, the equivalent param- 
eters are allowed to float, except no /z, or A c h parame- 
ters are included, and the K * (892) background yields are 
floated. 



A. LASS parameterization of (Ktt)q 

The J p = + component of the Kir spectrum, which 
we denote (Ktt)q, is poorly understood; we generate MC 
using the LASS parameterization 0, @ ; which consists of 
the A"q(1430) resonance together with an effective-range 
non-resonant component. The amplitude is given by 



cot 8b = 



rriKTT 



q cot 5b 

„2i8 B 



iq 



m r 



9o 



{ml - 
1 1 
aq 2 



K 



w )-im T - 



(8) 



(9) 



where man is the Kir invariant mass, q is the momen- 
tum of the Kir system, and go = l{ m o)- We use the 
following values for the scattering length and effective- 
range parameters: a = 2.07 ± 0.10 ( GeV) -1 and r = 
3.32 ± 0.34 (GeV) -1 For the resonance mass and 
width we use m = 1.412 GeV and r = 0.294 GeV. 

In the HMR, we parameterize the mKn distribution 
of the (-Ktt)o signal category with a Gaussian convolved 
with an exponential. This shape reasonably approxi- 
mates the LASS distribution, given the limited statis- 
tics in this analysis, and is chosen to reduce computation 
time. In the LMR, we use a linear polynomial, as only 
the tail of the (Kit)q enters the LMR wik-k region. 



VIII. BACKGROUND 5° 
YIELD FROM HIGH m,4 



/ 2 (1270)7r(892)° 
SIDEBAND 



To extract the f 2 (!270)K* (892)° yield, we select the 
LMR for mpcn and require a ir + ir~ invariant mass within 
the range 0.47 < m„+„- < 1.47 GeV. We perform an 
ML fit with the observables AE, tiies, and m^, 
and create a dataset of sWeighted if* (892)° events. We 
then fit the m- K +- K - spectrum of the sWeighted K*(892)° 
events to p°, f , and / 2 (1270) hypotheses (see Fig. [3]). 
We find 627±41 / 2 (1270)iT*(892) events after subtract- 
ing a 25 ± 13 event fit bias, which includes systematics; 
see Sec. HXI for details of the fit bias estimation method. 
The MC efficiency e of B Q -> / 2 (1270)if *(892)° is 11.8% 
(longitudinal polarization) and 20.4% (transverse polar- 
ization). 
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FIG. 3: (color online) -k + tv mass spectrum for sWeighted 
K*(892)° events. The solid curve is the fit function, the 
[green] dotted curve is p°, [red] dash-dotted is /o, and [blue] 
dashed is / 2 (127Q). 



B. PDF corrections from data calibration samples 

_ The decays B° -> D-ir+ (D~ -> K+ir-ir-) and B° -> 
D°ir° (D° — > K + ir~ir°) have the same particle content in 
the final state as the signal, as well as large branching 
fractions. They are used as calibration channels. We 
apply the same selection criteria described in Sees. Illlf ^ 
W\ except that the and rriK-n mass restrictions are 
replaced with 1.85 < m D - < 1.89 GeV or 1.83 < m^o < 
1.89 GeV and no D meson veto is applied. We use the 
selected data to verify that the ML fit performs correctly 
and that the MC properly simulates the J-, AE, and 
rriES distributions. From these studies, we extract small 
corrections to the MC distributions of AE and m^s, 
which we apply to the signal PDFs in our LMR and HMR 
likelihood fits. We find that it is not necessary to correct 
the PDF for T . 



Note that the three-component fit in Fig. [3] well de- 
scribes the resonances of interest, though the fit quality 
is poor at the lowest and highest 7r + 7r _ masses. As this 
study is intended to estimate the effect of higher reso- 
nances feeding into the nominal fit region, we determine 
that the three-component fit is sufficient. The excess of 
events in the low mass region could suggest the presence 
of a <r//o(600) resonance; this is accounted for in a sep- 
arate systematic study for the nominal fit. The excess 
in the highest bins could be explained by contributions 
from additional higher-mass resonances. As such reso- 
nances are unlikely to affect the p° and /o yields, we 
leave further understanding of these resonances for fu- 
ture studies. 

Using the MC efficiency for /^mO)^* (892)° in the 
LMR region, which includes a tighter cut on m n7r , 
and assuming fi = 0.5, we determine that there are 
47 ± 3 ^(1270)^* (892)° events expected in the LMR 
p°/f K* (892)° fit, as indicated in Table Q] 
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IX. FIT VALIDATION 



We include the isospin ratios 



Before applying the fitting procedure to the data, we 
subject it to several tests. Internal consistency is verified 
by performing fits to ensembles of simulated experiments. 
From these we establish the number of parameters asso- 
ciated with the qq PDF shapes that can be left free to 
float. Ensemble distributions of the fitted parameters 
verify that the generated values are reproduced with the 
expected resolution. 

We investigate possible biases on the fitted signal yield 
Yo, as well as on for the pK* (892) channels, due to 
neglecting correlations among discriminating variables in 
the PDFs, as well as from cross- feed from the BB back- 
ground modes. To determine these biases, we fit ensem- 
bles of experiments into which we embedd the expected 
number of signal and BB background events randomly 
extracted from detailed MC samples in which correla- 
tions are fully modeled. As correlations among fit vari- 
ables are negligible for qq events, these events are gener- 
ated from the PDFs. Each such experiment has the same 
number of signal and background candidates as the data. 
The measured biases are given in Table ITTll In calculating 
the branching fractions, we subtract the bias and include 
a systematic uncertainty (see Sec. IXIj) associated with 
the procedure. 

The two-stage fit employed to determine the 
p°/fo{Kn)f and /o^ 2 * (1430)° yields (see Sec. IVTBl) 
complicates the validation procedure. We perform the 
first stage of the fit (which extracts the "inclusive" 
(Ktt)q° and K^IHO) yields) on ensembles of experi- 
ments, as described above. The bias obtained from this 
study is split between the p° and fo channels based on 
the relative fraction of p° or fo events to the total number 
of signal events in that sample. 



X. FIT RESULTS 

The branching fraction for each decay chain is obtained 
from 



Y-Y 
eN B ft Bi 



(10) 



where Y is the yield of signal events from the fit, Yo is 
the fit bias discussed in Sec. lIXl e is the MC efficiency, Bi 
is the branching fraction for the i unstable B daughter 
(Bi having been set to unity in the MC simulation), and 
Nb is the number of produced B° mesons. The values of 
Bi are taken from Particle Data Group world averages @ • 
We assume the branching fractions of T(4S) to B + B~ 
and B a B° to be the same and to each equal 50%. As the 
branching fractions B(f — > itit) and B((Ktt)q — > Kit) 
are poorly known, we measure the products 

B{B° ->■ f K*) x B(f tttt) and 
B(B° -> X(Kir)*) x B((Ktt)* -> Kir) . 



r((Kn)f -> 


K+TT-) 


2 


T((Kn)f - 


-> Kir) 


3 


T((Kn)*+ 


> K+TT°) 
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r((Kn)* + - 


-> Kit) 


3 


r(/o(980) - 


> 7T + 7T~) 


2 


r(/ (980) 


— > 7T7r) 


3 



in our calculations of Jl^i- The efficiency e is eval- 
uated from the simulation. For the p~K* + channels, 
we apply an efficiency correction to the MC of roughly 
97%/7T°. The specific values are determined by calculat- 
ing a correction as a function of the ir° lab momentum 
from a detailed MC simulation of the signal channel. The 
correction is determined from a study of tau decays to 
modes with 7r°'s as well as a study of e + e~ — > 70; with 
ui — > 7r + 7T~7r°. The results for all signal channels are 
collected in Table Hill 

For all signals obtained from a one-stage ML fit, we 
determine the significance of observation S by taking the 
difference between the value of — 21n£ for the zero sig- 
nal hypothesis and the value at its minimum. For the 
P°(Ktt)^° : f (KTr)f, and f K^U3Q)° channels, the fit 
method does not readily provide a — 21n£ distribution, 
so we determine the significance assuming Gaussian un- 
certainties, which provides a conservative lower limit on 
S. 

For the p~(Ktt)q + channel, which has a significance 
less than 3ct including systematics, we quote a 90% C.L. 
upper limit, given by the solution Bgo to the equation 



j™C{b)db 



0.9. 



(11) 



where C(b) is the value of the likelihood for branching 
fraction b. Systematic uncertainties are taken into ac- 
count by convolving the likelihood with a Gaussian func- 
tion representing the systematic uncertainties. 

We show in Fig. [4] the data and fit functions projected 
onto the variable rriESi while in Fig. [5] we do the same for 
the 7T7T and Kir invariant masses for the LMR measure- 
ments. In Fig. EKa) we project the data and fit functions 
from the first stage of the HMR (Kn)^ Q and if 2 * (1430)° 
fits onto m K +T T -. In Figs. |H [SJ and Ufa) the signals are 
enhanced by the imposition of restrictions on the like- 
lihood ratio, which greatly reduce the amount of back- 
ground while retaining events that have a large probabil- 
ity to be signal. 

Figures [6jb) and (c) show the results of the second- 
stage HMR fit, distinguishing between the p° and fo hy- 
potheses. In these plots, we do not impose a restriction 
on the likelihood ratio, as these sWeighted samples al- 
ready contain only (b) (^vr)o or (c) ^2(1430)° signal 
events. 

Ref. Q extracts the resonant Xq(1430)° fraction of 
the LASS-parameterized (i^7r)Q° distribution. The reso- 
nant fraction is found to account for 81% of the LASS 
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TABLE III: Signal yield Y and its statistical uncertainty (see Sec. IVI 61 for an explanation of the two errors on the (Kty)*-, 
and 7f|(1430) yields); fit bias Yo; detection efficiency e for longitudinal (In) and transverse (tr) polarizations, if appropriate; 
daughter branching fraction product FJ Si ; significance S including systematic uncertainties; measured branching fraction B 
with statistical and systematic errors; 90% C.L. upper limit (U.L.); longitudinal fraction fi; and charge asymmetry A c h- In 
the case of foK*°, the quoted branching fraction is the product of B(B° — » foK*°) x B(fo — s* tytt). For the B° — > X(K-k)q 
channels, the quoted branching fraction is the product of B(B° — > X(Kty)q) x B((Kit)q —>■ Kty). We include the isospin ratios 
((Kir)? -> K+TT-) = 2/3, ({Ktt)*+ -> K + tt°) = 1/3, and (fo -> tt+tT) = 2/3. 



Mode 


y 

(events) 


Y 
(events) 


e(ln) e(tr) FJ^ 
(%) (%) (%) 


S 

(<r) 


B 

(io- 6 ) 


U.L. 

(IO" 6 ) 


h 




p°K* (892)° 
P°(Kny ° 


376 ± 37 
1045 ± 36 ± 118 


44 ±3 
80 ±11 


14.3 25.1 
9.6 


66.7 
66.7 


6.0 
6.3 


5.1 ± 0.6±g:| 
31 ± 4 ± 3 




0.40 ±0.08 ±0.11 


-0.06 ±0.09 ±0.02 


f K*(892)° 

h(K-K)l° 

/ A- 2 *(1430)° 


220 ± 23 
88 ± 19 ± 10 
134 ± 14 ± 23 


2.1 ± 1.6 
7± 1 
0±2 


18.3 
12.5 
15.3 


44.4 
44.4 
21.7 


9.8 
3.0 
4.3 


5.7 ±0.6 ±0.4 
3.1 ±0.8 ±0.7 
8.6 ± 1.7 ± 1.0 






±0.07 ±0.10 ±0.02 


p-A"*(892) + 
p-(JCir)S + 


167 ± 27 
221 ± 74 


23 ±3 
-5 ±8 


4.9 11.2 
4.5 


33.3 
33.3 


5.1 
2.8 


10.3 ±2.3 ±1.3 
32 ± 10 ± 6 


< 48 


0.38 ±0.13 ±0.03 


±0.21 ±0.15 ±0.02 




FIG. 4: (color online) _B-candidate iues projections for (a) 
p° I f K* (892)° (b) "inclusive" (Ktv)o° and 7^2(1430)°, (c) 
p-_R"*(892)+, (d) P -(K-k)* q + and p'KS (1430)+. The solid 
curve is the fit function, [black] long-dash-dotted is the total 
background, and the [blue] dashed curve is the total signal 
contribution. In (a) we separate the [red] dashed p° compo- 
nent from the [green] dotted fo- In (b) and (d) (Ktv)q signal is 
[green] dotted and 7^(1430) is [red] dashed. In (b), the two- 
stage nature of the fit means that the (Ktt)q° and K% (1430)° 
signals include both fo and p° components, as the first stage 
of the HMR fit does not include information about the tt + tt~ 
mass spectrum. 



shape in B + — > (Ktt)q°tt + decays. Using this reso- 
nant fraction along with the daughter branching fraction 
B(K£ (1430) -> Kit) = (93 ± 10)% @, we find the reso- 



FIG. 5: (color online) Invariant mass projections for LMR 
(a,c) p // jr(892) and (b,d) p~iY*(892) + ; tttv mass (left) 
and Kit mass (right). The solid curve is the fit function, 
[black] long-dash-dotted is the total background, and the 
[blue] dashed curve is the total signal contribution. In (a) 
we separate the [red] dashed p° component from the [green] 
dotted fo- In (c) and (d) (K-k)o background is [green] dotted. 



nant branching fractions 

B(B° -> p°K% (1430)°) = (27±4±2±3) x 10~ 6 ; 
B(B° -> /oAT *(1430) ) x B(f -> tttt) 

= (2.7 ±0.7 ±0.5 ±0.3) x 10~ 6 
B(B° -> p-iq(U30) + ) = (28±10±5±3) x 10" 6 , 

where the uncertainties are statistical, systematic, and 
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0.6 0.8 1.0 1.2 

nW(GeV) 

FIG. 6: (color online) Invariant mass projections for HMR 
p°(A"7r)5°, /o(A'7r)S°, and f K% (1430)° signals (a) K + tv~ 
mass, (b) ir + ir~ mass for sWeighted (Ktv)q° events, (c) tt + tt~ 
mass for s Weighted JsT|(1430)° events. The solid curve is the 
fit function. In (a) the [black] long-dash-dotted is the total 
background, the [blue] dashed curve is the total signal con- 
tribution, [green] dotted is the (Ktt)q° component, and the 
-R"2(1430)° is [red] dashed. In (b) and (c) the p° component 
is [red] dashed, fo is [blue] long-dashed. 



from the Xq(1430) — > Kit branching fraction, respec- 
tively. 



XI. SYSTEMATIC UNCERTAINTIES 

Tabic IIVI summarizes our estimates of the various 
sources of systematic uncertainty. We distinguish be- 
tween uncertainties that concern a bias on the yield (ad- 
ditive) and those that affect the efficiency and total num- 
ber of BB events (multiplicative), since only the former 
affect the significance of the results. The additive system- 
atic uncertainties are the dominant source of systematics 
for the results presented in this paper. The final row of 
the table provides the total systematic error in units of 
branching fraction for each channel. 



A. Additive systematic errors 

ML fit: We evaluate the systematic uncertainties due 
to the modeling of the signal PDFs by varying the rele- 
vant PDF parameters by uncertainties derived from the 
data control samples (see Sec. I VII B[) . This uncertainty 
is larger for the p~ K* + channels, as the D°7r° control 
sample has lower statistics than the D~ir + sample used 



for the p°K*° channels. 

Fit bias: The fit bias arises mostly from correlations 
among the fit variables, which are neglected in the ML 
fit. Studies of this bias are described in Sec. IIXI The 
associated uncertainty is the sum in quadrature of half 
the correction and its statistical uncertainty. For the 
p°/f (Kir)*° and f K* (1430)° channels, we add the un- 
certainty on the total bias in quadrature with half the 
bias scaled by the ratio of p° or fo events to their sum. 
BB background: We estimate the uncertainty from 
the fixed BB background component yields by repeat- 
ing the fit with the yields of these components varied 
by their uncertainties. For each signal channel, we add 
in quadrature the change in signal yield from varying 
each BB background, and quote this as the systematic. 
The uncertainty on the measured p(Ktt)q branching frac- 
tions makes this a large systematic for p°K* (892)° and 
p-K*(892)+. 

/o(980) parameters: The width of the fo is not accu- 
rately measured; to account for this, we allow the mean 
and width of the fo to float in the LMR fit and take half 
the shift in the signal yield as a systematic. This is one of 
the largest systematics for f K* (892)°. In the HMR, we 
lack the statistics to allow these parameters to float, and 
so perform the fit with them fixed to the parameters ob- 
tained when floating them in the LMR. This is amongst 
the largest systematics for / (X7t)q and p°(Ktt)q° . 
LASS shape parameters: For the (Ktt)q channels, we 
vary the LASS parameters in the MC by the uncertainties 
listed in Sec. I VII Al and re- fit the data sample with PDF 
parameters based on this new MC. In each channel, we 
take the largest variation in the yield as a result of this 
procedure as the systematic. The LASS systematic is the 
dominant one for p~ (Ktt)q + . 

p° / fo Interference: The interference between the p° 
and fo integrates to zero over the symmetric Hp range. 
Additionally, the differential rate is an odd function of 
H p , so the fact that we use \H P \ in the fit means that the 
interference term also vanishes from the differential rate. 
K* Interference: In our nominal fits, we do not ac- 
count for interference between the scalar and vector K* , 
or between the vector and tensor. We estimate the mag- 
nitude of the K* interference effect in a separate calcula- 
tion, which takes into account the relevant mass and he- 
licity acceptance functions, and varies the relative strong 
phases between components over the full range. As in- 
terference can affect the K* lineshape, we conservatively 
take this systematic to be additive. This is among the 
dominant systematic uncertainties in the HMR fits. 
/ 2 (1270) interference: The fit described in Sec. lVHTl 
used to estimate the background from B° — > 
/2(1270)-fr*(892)° decays, is performed without interfer- 
ence. Interference terms vanish when integrating over 
the full solid angle, however the requirement that the 
helicity angle be \H P \ < 0.9 leaves a non-zero interfer- 
ence term. For the case of interference between the fo 
and /2(1270), the scalar fo may interfere with the lon- 
gitudinal component of the /2(1270). Adding this term 
to the fit results shown in Fig. [3J assuming fa = 0.5 
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TABLE IV: Estimates of systematic uncertainties. 



Quantity 


v*0 

p K 


rs*0 


p ts*0 

foK u 


foK 


p TS-*0 

foK 2 


— 7^* + 

p K T 


P K o 


Additive errors (events) 
















ML fit 


2.7 


3.7 


1.1 


0.3 


0.7 


6.7 


21.3 


Fit bias 


ZZ.Z 


41. 


1.9 


Iz.o 


r, n 


1 1 n 
11.9 


Q 1 

0.4 


BB background 


14.8 


5.1 


2.7 


0.4 


0.6 


6.0 


3.2 


/o(980) parameters 


3.5 


10.0 


8.8 


10.0 


2.5 






LASS parameters 


— 


29.5 


— 


2.5 


7.7 


— 


31.0 


Interference 


12.2 


57.5 


12.9 


6.5 


9.9 


8.4 


18.1 


a// (600) 


+0.0 
-36.0 


+0.0 
-28.0 


1.0 


2.1 


6.4 


— 


— 


Total additive (events) 


+29.7 
-46.6 


+77.9 
-82.7 


16. U 


1 I .6 


1 c o 
lt>.2 


1 T 1 
1 1 .1 


42. 1 


Total additive [0(10 )J 


+0.45 
-0.71 


+2.58 
-2.74 


0.42 


0.67 


1.00 


1.24 


6.00 


Multiplicative errors (%) 
















Track multiplicity 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Track finding 


0.7 


0.7 


0.7 


0.7 


0.7 


0.4 


0.4 


Tr° efficiency 












3.8 


3.6 


Number BB 


0.6 


0.6 


0.6 


0.6 


0.6 


0.6 


0.6 


Branching fractions 










1.2 






IvlVy bLclLlbllCb 


u.uo 


U.UI 


n nfi 


U.U1 


n n^ 


n nfi 

u.uu 


n n^ 

u.uo 


COS 0T 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


1.5 


PID 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Jl uncertainty 


5.8 










2.1 




Total multiplicative (%) 


6.2 


2.3 


2.3 


2.3 


2.6 


4.9 


4.2 


Total systematic [B(10 _6v 


1 +0.6 
1 -0.8 


+2.7 
-2.8 


0.4 


0.7 


1.0 


1.3 


6.1 



for the /2(1270), and scanning over the unknown phase 
difference between the /o and /2(1270), we find a maxi- 
mum yield difference between the case of no interference 
of ±8.4% in the /o region. As the interference depends 
upon the sine of the unknown phase, we divide by y/2 
and report an additive / — / 2 (1270) interference system- 
atic of 12.8 events (5.8%). Using a similar procedure for 
p — /2(1270) interference, we report a systematic of 6.8 
events (1.8%). For the / iv*(892)° measurement, this is 
the dominant systematic. 

er//o(600) resonance: The scalar cr//o(600) is poorly 
understood and its parameters uncertain. We esti- 
mate the effect of a possible ct//o(600) resonance by in- 
cluding a/f (600)K* as a separate component in our 
fits. We parameterize the cr//o(600) using a relativis- 
tic Breit Wigner function with m = 513 MeV and T = 
335 MeV 0. As we lack a/f (600)K* MC, for the LMR 
K*(892)° fit we use the f K*° PDF shapes for all vari- 
ables except the 7r + 7r~ invariant mass. We use the aver- 
age cr// (600)iv*° branching fraction from the three K*° 
channels to calculate how many a/ fo(600)K*° events are 
expected in each K*° sample; this a/ f Q (600)K*° yield is 
then fixed in each fit. We take 100% of the resulting 
signal yield variation as a low-side systematic for the p° 
channels (a non-zero cr//o(600) yield decreases the p° 
yield) and conservatively consider this a two-sided sys- 
tematic in the /q channels. This is the dominant system- 
atic for p°K*(892)°. 



B. Multiplicative systematic errors 

Track multiplicity: The inefficiency of the selection 
requirements for the number of tracks in the event is a 
few percent. We estimate an uncertainty of 1% from the 
uncertainty in the low-multiplicity tail of the B decay 
model. 

Track finding/efficiency: Studies of tau events deter- 
mine that no efficiency correction is necessary for track 
finding and reconstruction. The systematic uncertainty 
is determined by adding linearly 0.17% per track in 
quadrature with an overall factor of 0.11%. 
7T° reconstruction efficiency: We apply an efficiency 
correction to the MC of roughly 97%/7r°; the correction 
depends on the ir° momentum spectrum, so is somewhat 
different in different channels. The uncertainty associ- 
ated with this correction is roughly 1.5%/ir . 
Number of BB events: A separate study [l8[ deter- 
mines the overall uncertainty on the number of produced 
BB pairs to be 0.6%. 

Branching fractions of decay chain daughters: 

This is taken as the uncertainty on the daughter particle 
branching fractions from Ref. Q. 

MC statistics: The uncertainty due to finite signal MC 
sample sizes (typically 430,000 generated events) is given 
in Table HI 

Event shape requirements: Uncertainties due to the 
cos 9t requirement are estimated from data control sam- 
ples to be 0.05 x (1 — (| cos6?t| cut value)). 



16 



PID: We estimate from independent samples that the av- 
erage efficiency uncertainty associated with particle iden- 
tification is 1.0%. 

fi, uncertainty: The signal yield reconstruction effi- 
ciency for VV channels depends on f^. As a result, any 
systematic uncertainty on ft, translates into a system- 
atic uncertainty on the efficiency through the following 
expression: 

- = f !V r M AfL (12) 

The systematic error on /x (A/z,) is given in Table fVl 

C. Charge asymmetry systematic errors 

From the analysis of a variety of data control samples, 
the bias on A c h is found to be negligible for pions and 
-0.01 for kaons, due to differences between K + and K~ 
interactions in the detector material. We correct the fit- 
ted A c h by +0.01 and assign a systematic uncertainty of 
0.02, mainly due to the bias correction. 

D. Systematic errors on /l 

Most systematic uncertainties cancel when calculat- 
ing /i. We include uncertainties from the signal PDF 
modeling ("ML fit"), fit bias (for which we assign an 
uncertainty equal to 100% of the bias added in quadra- 
ture with its uncertainty), BB background yields, the 
/o parameterization, and the possible existence of a 
cr//o(600) (where we take 100% of the Jl variation when 
the a/ fo(600)K*° is fixed in the study described in 
Sec. EH}. For p°K* (892)°, the fit bias of -0.045 ±0.008 
provides a moderate uncertainty; for p~K*(892) + , this 
bias is small (-0.009 ± 0.014). See Table|V]for details. 

TABLE V: Estimates of systematic errors on /l- 



Quantity 


p°iT (892)° p- 


-7T(892)+ 


ML fit 


0.003 


0.012 


Fit bias 


0.046 


0.016 


BB background 


0.019 


0.024 


/o(980) parameters 


0.004 




<t//o(600) 


0.100 




Total 


0.112 


0.031 



XII. DISCUSSION AND SUMMARY OF 
RESULTS 

We obtain the first observations of B° —> p (Ktt)q°, 
f Q K* (892)°, and p~K* (892)+ with greater than 5ct sig- 
nificance, including systematics. We present the first ev- 
idence for B° -> f Q (KTT)* ° and f K% (1430)°, which we 



observe with a significance of 3. Oct and 4.3ct, respectively. 
All branching fraction measurements have greater than 
3ct significance including systematics, except p~(Ktc)q + 
for which we also quote a 90% C.L. upper limit. No sig- 
nificant direct CP-violation is observed. Our results are 
consistent with and supersede those reported in Ref. @. 
For the if* (892) channels, we find the following results 



B(B° -> p°K* a ) 

h{p°K*°) 

A ch (p°K*°) 



p-K*+) 



B{B° 

A ch (p-K*+) 



(5.1 ± 0.6±g;|) x 10~ 6 
0.40 ±0.08 ±0.11 ; 
-0.06 ±0.09 ±0.02 



(10.3 ± 2.3 ± 1.3) x 10~ 6 
0.38 ±0.13 ±0.03 ; 
+ 0.21 ±0.15 ±0.02 ; 



B(B° -> f K*°) x B(/ -► tttt) 

= (5.7 ±0.6 ±0.4) x 10~ 6 ; 
AchifoK* ) = + 0.07 ± 0.10 ± 0.02 . 

The p°K*(892)° results agree with previous BaBAR @ 
and Belle [7J measurements and are consistent with pre- 
dictions from QCDF @. The p~K* (892)+ results are 
consistent with the previous BABAR upper limit and 
agree with QCDF predictions. Both the p° K* (892)° and 
p~K*(892) + branching fractions are, however, higher 
than the values predicted by QCDF. We find a branching 
fraction for J ^*(892)° within the previous BABAR 90% 
C.L. upper limit (6.5 x 10 -6 [6]) and somewhat above 
the Belle limit (3.3 x 10~ 6 7]), where we have scaled 
the published limits by a factor of 3/2, as the previ- 
ous analyses assumed B(/o — > 7r + 7r~) = 100% whereas 
this measurement includes the isospin ratio r(/o — > 
7T+7r-)/r(/ -> tttt) = 2/3. The / ^*(892)° branching 
fraction result is within one sigma of the QCDF pre- 
diction of 4.8^2,n x 10~ 6 , which is scaled by a factor 
of 3/4, as Ref. [| assumes B(f -> tt+tt-) = 0.5. We 
note that a previous BABAR study of B — > cj>K * Ql| ob- 
served an excess of B° -> (K + K~) K* (892)° events, 
where the scalar (K + K~)o could include / (980) de- 
cays. If we assume all the observed (K + K~)q excess 
to be from f -> K+K~ and follow Ref. @ in defining 
the ratio R = V (tttt) /[T (tttt) + T(KK)} ~ 0.75, then the 
B° -> f K* (892)° branching fractions are comparable 
for the /o — > tt + tt~ and K + K~ channels. 

As expected for penguin-dominated channels, the mea- 
sured Jl values are inconsistent with the naive factoriza- 
tion prediction of /l ~ 1. The predicted fz for p~K* + 
is higher than the measured value, though the theory er- 
rors are still large. Including the results from this paper 
and averaging BaBAR Q and Belle Q fh measurements 



*o 



we_can order the experimentally measured 
as 



for p+K 
values of ft 

h(p-K* + ) < f L (p°K*°) < f L (p + K*°) < f L (p°K* + ) , 

with the values ranging from 0.38 — 0.78. With the cur- 
rent experimental sensitivities, the three smallest Jl val- 
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ues are consistent with each other at la. QCDF [9] pre- 
dicts the following hierarchy among these fz values 

f L (p°K*°) < f L (p + K*°) < f L (p-K*+) < f L (p°K*+) , 

which agrees with the experimental finding that 
f l{p° K* + ) is largest. A more rigorous test of the theo- 
retical hierarchy requires additional experimental input. 
For fo(Kn)^ and /o-K3(143Q) , we find 

B(B° -+ .MKtt)* ) 

xB(/ -> tttt) x B{{Kn)l -> Ktt) 

= (3.1 ±0.8 ±0.7) x 10~ 6 ; 
B(B Q -> / tf 2 *(1430)°) x B(f -> tttt) 

= (8.6 ± 1.7 ± 1.0) x 10~ 6 . 

For p(Kit)q, we find, 

B{B° -> p a {KTr)*°) x B((#tt)S -> Xtt) 

= (31±4±3) x 10~ 6 ; 

B{B° -> p-(!f7r)* + ) x B{(Kn)* -> ^tt) 

= (32 ±10 ±6) x 10~ 6 ; 
< 48 x 10~ 6 . 

Using the ifo( 1430 )° resonant fraction of the LASS 
(Ktt)* result from Ref. ||, we can calculate the branch- 
ing fractions for the -K~q(1430) component of our (Ktt)q 
channels. We find 

B(B° ^ p°K*(U30)°) = (27±4±2±3) x 10~ 6 ; 
B(B° -> /o^ *(1430)°) x B(f -> tttt) 

= (2.7 ±0.7 ±0.5 ±0.3) x 10~ 6 
S(B° ->• p-^(1430)+) = (28 ± 10±5±3) x 10~ 6 , 



where the third uncertainty is from the daughter branch- 
ing fraction 6(^(1430) -> Ktt) = (93 ± 10)% 0. 
These results are somewhat lower than the QCDF pre- 
dictions @ but are consistent with QCDF within the un- 
certainties. The pQCD predictions have central values of 
(0.5 — 10) x 10~ 6 and are, in most cases, inconsistent with 
our results. 



XIII. ACKNOWLEDGMENTS 

We are grateful for the extraordinary contributions of 
our PEP-II colleagues in achieving the excellent luminos- 
ity and machine conditions that have made this work pos- 
sible. The success of this project also relies critically on 
the expertise and dedication of the computing organiza- 
tions that support BABAR. The collaborating institutions 
wish to thank SLAC for its support and the kind hospital- 
ity extended to them. This work is supported by the US 
Department of Energy and National Science Foundation, 
the Natural Sciences and Engineering Research Council 
(Canada), the Commissariat a l'Energie Atomique and 
Institut National de Physique Nucleaire et de Physique 
des Particules (France), the Bundesministerium fur Bil- 
dung und Forschung and Deutsche Forschungsgemein- 
schaft (Germany), the Istituto Nazionale di Fisica Nu- 
cleare (Italy), the Foundation for Fundamental Research 
on Matter (The Netherlands), the Research Council of 
Norway, the Ministry of Education and Science of the 
Russian Federation, Ministerio de Ciencia e Innovation 
(Spain), and the Science and Technology Facilities Coun- 
cil (United Kingdom) . Individuals have received support 
' from the Marie-Curie IEF program (European Union) 
and the A. P. Sloan Foundation (USA). 



[1] H.-Y. Cheng and J. G. Smith, Ann. Rev. Nucl. Part. Sci. 

59, 215 (2009). 
[2] Particle Data Group, K. Nakamura et al, J. Phys. G37, 

075021 (2010). 

[3] The mass and width of the fo (980) are from the Breit- 
Wigner parameters obtained by the E791 Collaboration, 
E. M. Aitala et al, Phys. Rev. Lett. 86, 765 (2001). 

[4] LASS Collaboration, D. Aston et al, Nucl. Phys. B 296, 
493 (1988). 

[5] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 72, 

072003 (2005). [Erratum- ibid. D 74, 099903 (2006)]. 
[6] BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

97, 201801 (2006). 
[7] Belle Collaboration, S. H. Kyeong et al, Phys. Rev. D 

80, 051103 (2009); Belle Collaboration, K. Abe et al, 

Phys. Rev. Lett. 95 141801 (2005). 
[8] BABAR Collaboration, P. del Amo Sanchez et al, Phys. 

Rev. D 83, 051101 (2011). 
[9] H. Y. Cheng and K. C. Yang, Phys. Rev. D 78, 094001 

(2008) [Erratum-ibid. D 79, 039903 (2009)]; M. Beneke, 

J. Rohrer and D. Yang, Nucl. Phys. B 774, 64 (2007); 

H. Y. Cheng, C. K. Chua and K. C. Yang, Phys. Rev. D 



77, 014034 (2008). 

[10] Z. Q. Zhang, Phys. Rev. D 82, 034036 (2010). 

[11] BABAR Collaboration, B. Aubert et al, Nucl. Instrum. 
Methods Phys. Res., Sect. A 479, 1 (2002). 

[12] The BABAR detector Monte Carlo simulation is based 
on GEANT4, S. Agostinelli et al, Nucl. Instrum. Meth- 
ods Phys. Res., Sect. A 506, 250 (2003) and EvtGen, 
D. J. Lange, Nucl. Instrum. Methods Phys. Res., Sect. 
A 462, 152 (2001). 

[13] S. Brandt et al, Phys. Lett. 12, 57 (1964); E. Farhi, 
Phys. Rev. Lett. 39, 1587 (1977). 

[14] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 79, 
072009 (2009). 

[15] M. Pivk and F. R. Le Diberder, Nucl. Instrum. Meth. A 

555, 356 (2005). 
[16] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. 

B 241, 278 (1990). 
[17] CLEO Collaboration, H. Muramatsu et al, Phys. Rev. 

Lett. 89, 251802 (2002). [Erratum-ibid. 90, 059901 

(2003)]. 

[18] G. McGregor, 0812:1954 [hep-ex] (2008). 

[19] BABAR Collaboration, B. Aubert et al, Phys. Rev. D 78, 



18 

092008 (2008). 



